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ABSTRACT 

We report the identification of elongated (triaxial or prolate) galaxies in cosmological 
simulations at z ~ 2. These are preferentially low-mass galaxies (M, ^ 10®'® Mq), 
residing in dark-matter (DM) haloes with strongly elongated inner parts, a common feature 
of high-redshift DM haloes in the ACDM cosmology. Feedback slows formation of stars 
at the centres of these halos, so that a dominant and prolate DM distribution gives rise to 
galaxies elongated along the DM major axis. As galaxies grow in stellar mass, stars dominate 
the total mass within the galaxy half-mass radius, making stars and DM rounder and more 
oblate. A large population of elongated galaxies produces a very asymmetric distribution of 
projected axis ratios, as observed in high-z galaxy surveys. This indicates that the majority 
of the galaxies at high redshifts are not discs or spheroids but rather galaxies with elongated 
morphologies. 
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1 INTRODUCTION 


The intrinsic, three dimensional (3D) shapes of today’s 
galaxies can be roughly described as discs or spheroids, or 
a combination of the two. These shapes are characterized by 
having no preferential long direction. Examples of galaxies 
elongated along a preferential direction (prolate or triaxial) 
are rare at z = 0 (iPadilla & Straussll2008t IWeiimans et al'l 
l2014l) . They are usually unrelaxed systems, such as ongoing 
mergers. However, at high redshifts, z = 1 — 4, we may 
witness the rise of the galaxy structures that we see today at 
the expense of other structures that may be more common 
during those early and violent times. 

Observations trying to constrain the intrinsic shapes of 
the stellar components of high-z galaxies are scarce but they 
agree that the distribution of projected axis ratios of high-z 
samples at 2 = 1.5 — 4 is inconsi stent with a population of 
randomly oriented disc galaxies (iRavindranath et al.l l2006t 


Law et al 


Law et al 


l2012t l^ma et alJl2012h . After some modeling, 
(12012l) concluded that the intrinsic shapes are 
strongly triaxial. This implies that a large population of 
hig h-z galaxies are elongated along a preferential direction. 

Ivan der Wei et al.l(l2014h looked at the mass and redshift 
dependence of the projected axis ratios using a large 
sample of star-forming galaxies at 0 < 2 < 2.5 

from CANDELS-I-3D-HST and SDSS. They found that the 
fraction of intrinsically elongated galaxies increases toward 
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higher redshifts and lower masses. They concluded that the 
majority of the star-forming galaxies with stellar masses 
of M* = 10® to 10®'® Mq are elongated at 2 ^ 1. At 
lower redshifts, galaxies with similar masses are mainly 
oblate, disc-like systems. It seems that most low-mass 
galaxies have not yet formed a regularly rotating stellar 
disc at 2 1. This introduces an interesting theoretical 

challenge. In principle, these galaxies are gas-rich and 
gas tends to settle in rotationally supported discs, if the 


angular momentum is conserved dFall & Efstathiou 

00 

iBlumenthal et al.ll986l;lMo et al.lll998l:lBullock et al. 

2001). 

At the same time, high-mass galaxies tend to be 

oblate 


systems even at high-z. The observations thus suggest that 
protogalaxies may develop an early prolate shape and then 
become oblate as they grow in mass. 


Prolateness or triaxiality are com mon properties of 


DM haloes in N-body-only simulations djing & Sut3 


All good et akl l2006t 


2008: Schneider et al 


Bett et al] 120071: iMaccio et al.l 


2001 


20071 


12012 , and references therein). Haloes 


at a given mass scale are more prolate at earlier times, 
and at a given redshift more massive halos are more 
elongated. Lor example, small haloes with virial masses 
around ~ 10^^ Mq at redshift 2 = 2 are as prolate as 
today’s galaxy clusters. Individual haloes are more prolate at 
earlier times, when haloes are fed by narrow DM hlaments, 
inc luding mergers, rather than isotropically, as described 
in IVerakliro et al.l (1201 111 . The progenitors of MW-size 
haloes are fairly prolate at redshift 2 = 2 and they are 
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increa singly more elongated at smaller radii (lAllgood et al.l 
I 2 OO 6 I 1 because their inner shells collapsed earlier. 

The shape of the inner D M halo could influence th e 
shape of the central galaxy (iDekel & Shlosma^ Il983h . 
If a triaxial halo dominates the inner gravitational 
potential, the inner galaxy feels a non-axisymmetric, 
non-rotating potential. Under these conditions, stellar 
orbits elongated along the major axis of the triaxial 
halo are c ommon within some char acteristic radius of 
the system (iBinnev & Tremain^I2008L §3.3.1 and 3.4.1). 
Therefore, elongated galaxies are expected within the 
ACDM paradigm. 

The fact that more massive dark matter haloes 
tend to be more prolate should be reconciled with the 
observation that more massive galaxies, presumably living 
in mo re massive haloes, are more disky (Ivan der Wei et akl 
12014 . The missing piece in the picture of elongated 
galaxies is the counter-effect of baryons on the shape 
of the inner halo and its corresponding gravitational 
potential. N-body-i-gasdynamics simulations of galaxy 
forma tion could proyide this missing link. IZemp et aU 
(I 2 OI 2 I) compared the shape of z = 2 dark matter haloes 
in adaptive mesh refinement (AMR) simulations with and 
without cooling, star formation and feedback. The addition 
of galaxy formation processes led to rounder haloes, 
especially in their inner parts, where baryons dominate 
the total mass. However, these simulations had an excess 
of baryon s compared t o estimates based on abun dance 
matching (iMoster et al.l 120131: iBehroozi et al.l 12013l) such 
that the simulated galaxies were always oblate, and the 
effect of baryons on the dark-matter halo was possibly 
overestimated. 

The Eris sim ulation of a galaxy with = 7.9 x 
10^^ Mq at z = 0 (iGuedes et al]l201ll) shows more realistic 
baryonic fractions at z = 0 and it shows the formation of 
a pseudobulge at z = 4 triggered by minor mergers tha t 
induced non-axisymmetric instabilities (iGuedes et al.l201^ . 
The mass in this early pseudobulge was about IO^Mq, 
within the mass range of elongated galaxies. Indeed, its 
stellar surface density map at z = 4 looks marginally 
triaxial. However, there is no inform ation about the DM 
mass and shapes in the inner regions. iButskv et al.l (12015h 
find elongated halo shapes in N-body-only and hydro 
simulations at z = 0 when the virial mass M^ir ^ IO^^Mq, 
while for higher M^ir the hydro simulations show 
substantially rounder halo shapes. However, they do not 
discuss the stellar shapes of their simulated galaxies or 
extend their analysis to higher redshifts, when haloes 
with low stellar-to-h alo mass ratios are more elongated 
(lAllgood et al.ll200^ . Thus, previous simulations have not 
revealed elongated galaxies as observed at z ^ 1. 

In this paper, we report the formation of elongated 
galaxies in cosmological hydrodynamic simulations. We aim 
at understanding why elongated galaxies form preferentially 
in low-mass haloes at high redshifts. We also address 
whether the gas distribution is also elongated and whether 
the simulations can reproduce the typical shapes observed 
in the distribution of projected axis ratios in high-z galaxy 
surveys. In section Sj2]we describe the sample of simulations. 

describes the formation of elongated galaxies within 
prolate haloes and @ discusses their transformation into 
disk or oblate galaxies due to a high stellar central density. 


Table 1. Properties of the zoom-in simulations at 2 = 2. Columns show 
the name of the run, the virial radius, the virial mass, the stellar half-mass 
radius, the stellar mass and gas mass. Masses are in Mq and radii in kpc. 


Run 

Rv 

Mv 

ro.s 

Stellar mass 

gas mass 

Vela05 

44 

0.7 x toll 

1.8 

0.7 x 10® 

0.5 x 10® 

Vela04 

53 

1.2 x lO^i 

1.7 

0.8 x 10® 

0.8 x 10® 

Vela02 

55 

1.3 X loll 

1.9 

1.6 x 10® 

1.2 x 10® 

Vela28 

63 

2.0 X lOii 

2.3 

1.8 x 10® 

2.1 x 10® 

VelaOl 

58 

1.6 X lOii 

0.9 

2.0 x 10® 

1.5 x 10® 


SjS] shows the distribution of projected axis ratios and Sj6] 
summarizes the results and discusses them. 


2 SIMULATIONS 


The zoom-in simulations used in this paper were drawn from 
a larger dataset, the Vela sample, which contains 35 haloes 
with masses between 2 x 10^^ Mq and 2 x 10^^ Mq at 
z = 1. We selected the haloes that contain a galaxy with 
a stellar mass lower than 10®'^ Mq at z = 2. These are 
the best candidat es for hosting elongated galaxies, according 
to observations (Ivan der Wei et al.ll20l4 . The properties at 
z = 2 of the five haloes that meet this criterion are described 
in Table[T] Their typical virial mass is ~ 10^^ Mq and 
each hosts a galaxy with a stellar mass of M* ~ 10® Mq, 
0.6 — 1.2% of the total virial mass. This M*/Mv ratio is 
within a factor 2-3 of the estimates with abundance matching 
techniques, which are highly unc ertain at this redshift and 
mass range (IBehroozi et akluM^ . 

The simulations were performed w ith the ART code 
(iKravtsov et al.l 119971: iKravtsovI l2003ll . which accurately 
follows the eyolution of a grayitating N-body system and 
the Eulerian gas dynamics using an AMR approach. Beyond 
gravity and hydrodynamics, the code incorporates many 
of the physical processes relevant for galaxy formation. 
These processes, representing subgrid physics, include 
gas cooling by atomic hydrogen and helium, metal and 
molecular hydrogen cooling, photoionization heating by 
a constant cosmological UV background with partial 
sel f-shielding, star formation and feedback, a s desc ribed 


ICeverino & Klvpin (l2009l) . ICeverino et akl (1201 Ol) . and 
ICeverino et al.l (12014 ). In addition to thermal-energy 


feedback, the simulations discussed here use radiative 
feedback. This model adds a non-thermal pressure, radiation 
pressure, to the total gas pressure in regions where ionizing 
photons from massive stars are produced a nd trapped. In 
the cur rent implementation, named RadPre in ICeverino et al.1 
(I20l4 . radiation pressure is included in the cells (and their 
closest neighbors) that contain stellar particles younger than 
5 Myr and whose gas column density exceeds 10®^ cm“®. 

The initial conditions of these runs contain 6.4 — 11 x 
10® dark matter particles with a minimum mass of 8.3 x 
10^ Mq, while the particles representing stars that were 
formed in the simulation have a minimum mass of 10® Mq. 
The maximum spatial resolutio n is between 17-35 pr oper 
pc. More details can be found in ICeverino et al.l (12014 and 
IZolotov et al.l (l2015l) . 
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Figure 1. Example of an elongated galaxy. Orthogonal views along the galaxy principal axes, showing gas, stars, and DM (from top to bottom) of Vela28 
at 2 : = 2.2. The size of the images is 20 kpc. The contours of stellar surface density are defined as log{E/(M 0 kpc~^) = 1, 1.5 in the b-a and c-a planes 
and as log(E/(MQ kpc“^) = 1.5, 2 in the b-c plane. The stellar half-mass radius is ro .5 = 2.1kpc. The contours of DM surface density are set as 
log(E/(MQ kpc~^) = 2, 2.5 along the minor and intermediate axes and as log{E/(MQ kpc“^) = 2.2, 2.8 along the major axis. The labels at the top-right 
corner mark the directions of the galaxy axes in each view. 


3 FORMATION OF ELONGATED GALAXIES 

Figure [T] shows three orthogonal views of the three 
components (gas, stars, and DM) of an elongated galaxy at 
z = 2.2 (Vela28). The face-on view (b-a plane) was selected 
using the angular momentum of the cold (T < 10^ K) gas 
inside 5 kpc. The gas shows a very irregular and multiphase 
medium with clouds of relatively high column densities, 
log(I]/M 0 pc“^) = 1.5 — 1.7, and bubbles with very low 
densities. The geometry is far from a uniform, thin disc. 


although the gas is flattened along a minor axis, c, coincident 
with the rotation axis of the gas. Indeed, the gas rotates at 
Vrot = 80 kms“^, averaged over a 2-5 kpc shell. The radial 
velocity dispersion measured at 200 pc scales is also high, 
(Tr = 30kms ~^, typical of very pe rturbed, turbulent and 
thick discs. See lCeverino et al.] (12012h for more details. 

The stellar distribution is very different from the thick 
gaseous disc. It is elongated along one direction, which 
defines the galaxy major axis a. The 3D isodensity surface 
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crossing the major axis around the half-mass radius is well 
fitted by an ellipsoid with axial lengths equal to (a, b, c) = 
(2.2, 0.76, 0.67) kpc. The intermediate-to-major axis ratio is 
low, b/a = 0.35, inconsistent with an axisymmetric disc. 
The shape is close to a prolate ellipsoid. 

The dark matter also shows a prolate shape with 
similar properties at roughly the same radius, (a, b, c)dm = 
(1.8, 0.73, 0.6) kpc and (6/a)DM = 0.4. The fact that the 
dark matter is also elongated in the same direction as the 
stellar component is relevant to the formation of such an 
elongated galaxy. 

For a better understanding of the connection between 
the stellar and DM shapes, we perform a quantitative 
analysis of the intrinsic shapes of the stellar and DM 
distributions. We select the major progenitors in snapshots 
between z = 7 and 1. At each snapshot, we determine the 
shapes by fitting ellipsoids to sets of particles ha ving (nearly) 
the s ame nearest neighbours-based density (iJing & Sutol 
12003) . we first compute the stellar density at the position of 
each stellar particle using a 80-particles kernel. The density 
at a given point is computed as the mass of the 80 closest 
particles divided by the volume of a sphere with a radius 
equal to the distance to the 80th nearest neighbour. 

We compute a first guess of the major axis direction 
at a radius r in the following way. We define a spherical 
shell at that radius with a thickness of 6r = 0.01 kpc. We 
compute the three eigenvectors of the inertia tensor of all 
stars within that shell. The eigenvector corresponding to the 
smallest eigenvalue defines the first guess of the direction of 
the major axis. 

We define a 3D isodensity surface that intersects the 
direction of the major axis at a radius r. The surface is 
selected by a volumetric stellar density (ps), defined as the 
average of the densities of the two points, located along each 
side of the major-axis at the intersection with the sphere of 
radius r. We also compute its standard deviation (as). The 
isodensity surface is traced by all particles with densities 
around the selected value, {p = Ps ± CTs)- In most cases, 
iTs < Ps and the density declines so fast with radius that the 
thickness of the surface is small, and we find that it is well 
fitted by a triaxial ellipsoidal shell. 

Finally, an ellipsoidal fit to that isodensity surface 
defines the 3 main axes of the ellipsoidal shape. The fit is 
done by diagonalizing the inertia tensor (normalized to unity 
mass) defined by all particles tracing the isodensity surface. 
The eigenvalues (ei < 6 2 < e.Q defin e the 3 main axes of a 
triaxial ellipsoidal shell (lRouthll2013L chapter 1) as 

a = (1.5(63-I-62 - ei))°-^ (1) 

b = (1.5(63 + 61-62))°'® (2) 

c = ( 1 . 5(62 + 61 - 63 ))°'® (3) 

where a is the major axis, b is the intermediate axis, and 
c is the minor axis of the ellipsoid. A new 3D isodensity 
surface is defined along the direction of the current major 
axis at a radius r. The process repeats until convergence is 
reached, usually after 2-3 iterations, when the values of the 
principal axes do not change significantly. The final results 
are profiles of the stellar shape at different radii. A global 
shape is defined as the average over all shells around the 
stellar half-mass radius, r = ro.5 ± 0.1 tq.s, where the 
half-mass radius, r-Q.s, is defined using all stars within 0. Irv. 


The same procedure is used for dark matter. The average 
shape of the DM halo is taken at the same stellar half-mass 
radius. Therefore, the intrinsic shape of stars and dark 
matter are computed at the same galac tocentric radius. The 
triaxiality parameter (iFranx et al.lll99lh is defined as: 



(4) 


The usual convention defines a shape as oblate if 0 ^ 
Tri < 1/3. A prolate shape is defined by a high value, 
2/3 < Tri ^ 1, while triaxial shapes have intermediate 
values, and triaxiality is not defined for a sphere. 

The top-left panel of Figure |2] shows a one-to-one 
equivalence between the elongation parameter for stars and 
dark matter: 

6 = y/T^^JbJ^. (5) 

This correspondence confirms the impression of Figure [T] 
At the galaxy half-mass radius, stars and dark matter share 
the same intermediate-to-major axis ratio. We find that if the 
DM shape is strongly elongated, high e and low b/a, stars 
are also elongated in the same direction (as in Figure [U. 
At 3 ^2^1, the simulated galaxies of low masses, 
log(M*/MQ) ^ 9.5, show a large fraction of prolate 
galaxies, Tri > 2/3, with a high triaxiality value of Tri = 
0.8 ± 0.2, a high elongation parameter, e = 0.8 ± 0.1, 
and a relatively low intermediate-to-major axis ratio, b/a = 
0.6 ± 0.2. In these low-mass, prolate galaxies, dark matter 
dominates the mass (and the gravitational potential) inside 
the galaxy half-mass radius. The median ratio between stars 
and dark matter masses inside tq.s is 0.3 ± 0.3. Thus, 
the formation of an inner elongated DM halo in the early 
Universe, as predicted by ACDM, gives rise to galaxies 
elongated in the same direction. 


4 THE END OF ELONGATION 

Figure |2] also shows the transition from low-mass prolate 
galaxies to more massive less elongated galaxies at lower 
redshifts, 1 ^ z ^ 2. Galaxies with M* > 10° ® Mq 
have a much lower elongation parameter, e = 0.4 ± 0.2, 
which corresponds to a high axis ratio, b/a = 0.9 ± 
0.1. Therefore, these more massive galaxies are close to 
axisymmetry around the minor axis. 

The key variable controlling this lack of symmetry 
seems to be the stellar-over-DM mass ratio inside ro.5. 
The gas mass within rg.s is lower or similar to the stellar 
mass. The median axis ratio for dark-matter dominated 
galaxies is low, b/a = 0.6 ± 0.2, similar to the values of 
low-mass and prolate galaxies. The galaxies dominated by 
stars, M,(ro,5 ) > Mdm(»"o. 5 )j have 6/a = 0.9 + 0.1, similar 
to the values of higher-mass, near axisymmetric galaxies. 
As galaxies grow in mass, the baryons concentrate at the 
centres of their halos, making stars and dark ma tter rounder 
and m ore axisymmetric around the minor axis (IZemp et al.l 
l2012h . which usually coin cides with the direction of the 
galaxy angular momentum (ICeverino et al.ll2M3120151 and 
references therein). 
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Figure 2. Con'elations between the elongation parameter, e, for star's and dai'k matter, the triaxiality pai'ameter, Tri, and the ratio of stars and DM masses at 
ro.5- There is a one-to-one equivalence between the elongation parameter for stars and dark matter, so that there is a large population of elongated, prolate 
galaxies (Tri > 2/3) especially among low-mass galaxies, M, ^ 10®'® Mq, dominated by dai'k matter M*(ro, 5 ) < MjjiyrCfo.s)- 


5 PROJECTED AXIS RATIOS 

Our current sample of simulations is too small for a 
direct and quantitative comparison with the distribution of 
projected axis ratios in current galaxy surveys. However, 
we can make a qualitative comparison and check whether 
elongated galaxies can generate asymmetric distributions 
of projected axis ratios, as observed in current samples 
of low-mass galaxies at redshifts 0 < z < 2.5 

(Ivan der Wei et al.]l2014l) . 

In this exercise, we take one of the simulations in our 


sample, Vela04, which shows an elongated morphology, 
similar to the example in Figure [T] Between z = 2.5 and 
2, it has an elongation parameter of e ~ 0.9 {b/a = 
0.4). At later times, the galaxy is less elongated, e ~ 0.6 
(h/a = 0.8) because stars dominate the mass within the 
half-mass radius, M,(ro.5)/MDM(To.5) — 2, as opposed 
to M,(ro,5)/MDM(To.5) — 0.4 at earlier times (2.5 < 
z < 2). We select 10 snapshots between z = 2.5 and 
z = 1.5 and measure the projected axis ratio in 1000 
random lines-of-sight per snapshot. The projected axis ratio 
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0.49 1.01 

q (projected axis ratio) 

Figure 3. Distribution of projected axis ratio of Vela04, using 10000 
different lines-of-sight from 10 snapshots between z = 2.5 and z = 1.5. 
The asymmetric distribution, biased towards low axis ratios, is typical of 
intrinsically elongated galaxies. 


is measured by fitting a 2D isocontour of surface density at 
ro.5 with an ellipse of axis ratio q. 

The distribution of projected axis ratios (Figure O is 
clearly asymmetric. It is biased towards low projected axis 
ratios, with an asymmetric peak around q ~ 0.35, a sharp 
drop below q = 0.3, and a long tail towards g = 1. 
This shape is qualitatively similar to the distribution shown 
in 3D-HST+CANDELS data for \o^(MJMcr,) = 9 - 1 0 
at 1.5 < z < 2.0 (Figure 1 in Ivan der Wei et ^ l2014t) . 
This comparison supports the idea that the asymmetry in 
the distribution of projected axis ratios is an indication of 
a large population of intrinsically elongated galaxies among 
low-mass galaxies at high redshifts. 


6 DISCUSSION AND CONCLUSIONS 

We addressed the formation of elongated galaxies in the 
early Universe using five zoom-in AMR cosmological 
simulations of low-mass galaxies, 10® ® Mq ^ ^ 

10® ® Mq at z = 2. Elongated galaxies form preferentially 
in highly elongated, low-mass haloes at high redshifts. 
This is due to asymmetric accretion from narrow filaments 
and mergers. At the same time, feedback prevents the 
overabundance of baryons at the centre of these small 
haloes. As a cautionary note, different feedback models 
may produce different galaxy shapes. Eor example, 
simulations witho ut radiative feedback overproduce stars 
by a larger factor (ICeverino et al.ll2014t iMoodv et al.ll2014h 
and generate much less prolate galaxies, according to their 
H-band mock images (Tang, Primack, et al. in preparation) 
and 3D shapes (Tomassetti, Sai, et al. in preparation). 
The gravitational potential is dominated by the elongated 
dark matter, so that the potential is highly prolate. This is 
consistent with prolate (Tri ~ 0.8) galaxies, elongated in the 
direction of the DM major axis (e = 0.8 or b/a ~ 0.6). 

As galaxies grow in mass, baryons start to dominate 
the central potential. This makes the inner halo rounder, 
probably due to the deflection of elongated orbits by a 


central mass concentration (e.g. iDebattista et al.ll2008h . We 
found a decline in the fraction of elongated galaxies when 
stars dominate the mass inside the half-mass radius. These 
galaxies are near axisymmetry and oblate (b/a ~ 0.9). 
These massive galaxies build a disc-Hspheroid embedded in 
a ro under halo, as descri bed in ICeverino et al.l (1201 Ol 120151) 
an d IZolotov et al.l (l2015l) . __ 

The process of wet compaction (iDekel & Burkertl[2014l: 
IZolotov et al.ll2015l) . in which a gas-rich galaxy shrinks into 
a compact star-forming spheroid, could drive the transition 
from a DM-dominated elongated galaxy to a baryonic 
self-gravitating spheroid. Mergers, counter-rotating streams, 
and the associated violent disc instability drive baryons to 
the galaxy centre, increasing the central stellar-to-DM mass 
ratio, making DM and baryons rounder. The transition from 
DM to baryon dominance ro ughly occurs at the po int of peak 
compaction (Eigs. 2 and 3 of IZolotov et al.llToiSl) . 

Although the stellar distribution is elongated, gas is 
distributed in a turbulent and thick rotating disk. These 
different shapes may be related to the different angular 
momentum content o f gas and stars . Acc ording to the 
framework outlined in iDanovich et al.l (l2015l) . accreted gas 
has a high spin, mostly coming from angular momentum 
acquired in the cosmic web. In the inner halo, gas loses 
some of this angular momentum due to internal torques. 
Moreover, if the inner halo is strongly elongated, the transfer 
of angular momentum from baryons to dark matter is intense 
such that stars may lose almost all their spin and rotational 
support, because they spend more time in the inner halo 
than the freshly accreted gas. This could explain why gas 
forms a rotationally supported disc, while stars have an 
elongated, non-rotating shapes. At later times, in more 
massive galaxies, their inner haloes are rounder, minimizing 
the angul ar momentum transfer and en abling stellar discs, a s 
shown in ICeverino et al.l (12015l) and in ISnvder et akl (l2014l) . 

Intrinsically elongated galaxies may generate a very 
asymmetric distribution of project ed axis ratios, as observe d 
in low-mass galaxies at z ^ 1 (Ivan der Wei et al.l l2014t) . 
Lower-mass galaxies are always more abundant than 
high-mass galaxies at any redshift. This implies that the 
majority of the galaxies at these redshifts are not discs 
or spheroids, but low-mass elongated galaxies. If such 
galaxies c ontribute decisiv ely to the reionization of the 
Universe (llliev et akl 12012l and references therein), this 
elongated morphology may have important implications for 
reionization. The ability of ionizing photons to escape their 
host galaxies may depend on the actual galaxy geometry. 
The geometry of galactic outflows may also depend on 
the degree of elongation of its host galaxy, because wind 
material tend to escape along the direction with less 
resistance. Therefore, very elongated galaxies may host 
outflows with cylindrical or toroidal symmetry, as opposed 
to the conical outflows expected in large discs. 
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